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Abstract 

Sterile neutrinos are one candidate to explain anomalies in neutrino oscillations. The mass- 
difference-driving oscillation between flavors can be probed only within specific combinations of 
baseline and flight energy. For a neutrino whose mass is completely unknown, it is necessary to scan 
all available ranges in spectrum and all accessible baselines. Here, we present four-neutrino analysis 
of the results announced by RENO and Daya Bay, which performed the definitive measurements of 
013 based on the disappearance of the inverse-beta-decay antineutrino at km-order baselines. Our 
results within 3+1 scheme include the exclusion curve of Am4i vs. On, and the adjustment of 6*13 
due to the contribution of ^14 to the disappearance of electron antineutrinos. 
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I. INTRODUCTION 



On top of the solar and atmospheric neutrino experiments, a series of recent oscil- 
lation experiments, T2k[i|, MINOsQ, Double ChoozS^, Daya Bay|4|, and RENofl, 
have firmly established a framework of neutrino oscillations among three flavor neutri- 
nos mixed with three mass eigenstates through a unitary matrix. According to the Su 
global analysisjo], the 3a ranges of the physical oscillation parameters are given as in 
the following: 6.99 < AmlJlQ-^eV^ < 8.18, 2.59 < sin^^is/lO-^ < 3.59, 2.19(2.17) < 
Amis/lO-W^ < 2.62(2.61), 1.69(1.71) < sin^^ig/lQ-^ < 3.13(3.15), and 3.31(3.35) < 
sin^ 6'23/10^^ < 6.36(6.63) for normal (inverted) hierarchy. While all the three mixing angles 
are now known to be different from zero, the values of CP violating phases are completely 
unknown yet. Although there are a number of global analysis which presented the values of 
masses and mixing parameters consistent with among themselves 6|-[8|, we focus on and 
its relatives obtained by RENO and Daya Bay. 

In spite of the confirmation of three flavors of neutrinos, we do not concretely exclude the 
existence of new kinds of neutrinos. According to the LEP experimental result for invisible Z 

n 

boson decay [9l], if there exist new types of neutrinos with mass below 45 GeV, they should 
be sterile neutrinos which are singlet fields under the SU(2) weak interaction. Although 
sterile neutrinos do not interact with the electroweak gauge bosons, they can mix with three 
active neutrinos, leading to the oscillation between active and sterile neutrinos. 



Inactive singlet neutrinos are familiar in utilizing the see-saw mechanism lOj, but they 



are untouchably heavy. On the other hand, the existence of light sterile neutrinos with 
masses about 0(1) eV or less has not been phenomenologically ruled out. It is naturally 
considered that their existence may affect cosmology such as Big Bang Nuclearosynthesis, 
Cosmic Microwave Background, Hubble constant and galaxy power spectrum etc. Many 
literatures have studied the impact of sterile neutrinos on cosmology and obtained some 
constraints on the effective numbers of light neutrino species and on the sum of light neutrino 
masses which may favor the existence of sterile neutrinos lll|. The anomalies observed in 



the LSND MJ, MiniBooNE 



13| . Gallium solar neutrino experiments [l^ and some reactor 



experiments 15| , over the past several years, also can partly be reconciled by the oscillations 
between active and sterile neutrinos, if more than one kind of sterile neutrino are heavier 
tha. t.n.ee active neutn.os H . 
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We examine whether the oscillation between a sterile neutrino and active neutrinos is 
plausible, especially by interpreting the results released from Daya Bay and RENO. For the 
sake of simplicity, we assume that one type of sterile neutrino is added into the contents of 
neutrinos. The analysis is restricted within a narrow range of Am^j^, since the setups of two 
experiments are optimized for Amg^. The search for Aml^ oscillation with a sterile neutrino 
is possible only if the order of Amli is not far different from the order of Am|^. There 
is no hope to probe a sterile neutrino in kilometer-baseline reactor neutrino oscillations 
such as RENO and Daya Bay, if the mass squared difference between the sterile and active 
neutrinos is not less than leV^. Thus, our study does not necessarily cover the sterile 
neutrinos introduced to reconcile the anomalies observed in the short baseline experiments. 
In fact, efforts to search sterile neutrinos are being made with all types of oscillations with 
different baselines 17H25|. There is a work that tried to probe sterile neutrino parameters 



with recent reactor neutrino experiments [20|, |21[ , where a different range of Ami^ is focused 
on, compared with this work. 

This article is organized in the following outline: In Section II, the survival probability 
of electron antineutrinos is presented in four- neutrino oscillation scheme. We exhibit the 
dependence of the oscillating aspects on the order of Amli, when the reactor neutrinos with 
the energy range 1 to 8 MeV are detected after a travel along km-order baseline. In Section 
III, the rate-only analysis results announced by RENO and Daya Bay are re-analyzed in 
version of four-neutrino oscillation, and the relevant spectral shape analysis follows in the 
next section. In Section IV, the curves of four-neutrino oscillation are compared with the 
data obtained at the experiments in order to search for any clue for a sterile neutrino and 
in order to see the change in sin^ 26'i3 in the coexistence with the sterile neutrinos. Broad 
ranges of Amli and sin^ 29i4 remain not being excluded. The exclusion bounds of sin^ 29i4 
and the best fit of sin^ 2^13 are summarized in conclusion. 

II. FOUR NEUTRINO OSCILLATION 

The three-neutrino transformation from mass basis to flavor basis is given in terms of 

I — I 

three angles and a Dirac phase [29,] : 

^PMNS = -R23(^23)-Rl3(^13, '^l)-Rl2(^12), (1) 
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where RijiOij) denotes the rotation of the i-j block by the angle 6ij. When a 3+1 model 
is assumed as the minimal extension, the unitary transformation from the mass basis of 
{mi, m2, m3, 7/14} to the flavor basis {i^e, ^^t, ^s} is given in terms of six angles and 
three Dirac phases: 

Uy = -R34(6'34)-R24(6'24, ^2)-Rl4(^14) • 

■R23{023)Rl3{0l3,Si)Ru{eu,63). (2) 

The 4-by-4 Up is expressed as 



I f/el f/e2 f/e3 0\ 

U^l U^2 U^3 

Url Ur2 Ur3 

\ 1/ 



(3) 



(4) 



/ Ci4 Si4 \ 

— S14S24 C24 C14S24 

— C24S14S34 — S24'534 C34 C14C24S34 
V-C24C34S14 -S24C34 -S34 C14C24C34/ 

^CiiUel Ci4t/e2 Cl4Ue3 "5l4 ^ 

C14S24 

C14C24S34 

V C14C24C34/ 

where the PMNS type of 3-by-3 matrix ?7pmns with three rows, (f/ei f^e2 Ue3), {U^i U^2 U^^^) 
and {Uri Ur2 Ur3), Is imbedded. The CP phases 62 and ^3 introduced in Eq.([2]) are omitted 
for simplicity, since they do not affect the electron antineutrino survival probability at the 
reactor neutrino oscillation. 

The survival probability of Ue produced from inverse beta decay is 



exp t- 



\u, 



i<j 



(5) 
(6) 



where Amfj denotes the mass-squared difference rnf ""^j- the limit where \Ue3\ and \Ue4\ 
are much smaller than one, the survival probability of Ue is determined by the following 
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three terms (as long as is much larger than the others): 



Pth{i^e i^e) = 1 — cf4Ci3 sin^ 29i2 sin^(1.27Am: 




- sin^ 2^13 sin2(1.27Am^^—) 



E 



(7) 



— sin 2^14 sin2(1.27Am2^— ). 



E 



The oscillation pattern of Pth as L/E increases in a logarithmic way is described in Fig HI 
where three patterns of oscillating probabilities are shown according to the order of Am4^. 
It is shown that the order of /S.m\i to probe a 4th neutrino must be not much larger than 
that of Amg^, since both RENO and Daya Bay have baselines optimized for Am\^. The 
first bump in each curve corresponds to the oscillation due to Am\^, while the second bump 
that appears near 4 ~ 5m/MeV corresponds to the oscillation due to Am|]^. If ^m^i 
is much less than (9(0.01)eV^, the amplitude of the Amli oscillation can appear only in 
the superposition with the Amg^^ oscillation. The bound on Aml-^ at which the two types 
of bumps can separate is about O.OOSeV^. On the other hand, if Amli > (9(0.03), the 
oscillating aspect of the far-to-near ratio becomes unclear. Hereafter, we consider the mass 
of the sterile neutrino m4 within O.OOSeV^ < Amj^ < O.OSeV^ 

III. RENO AND DAYA BAY: REACTOR NEUTRINO OSCILLATION EXPERI- 
MENTS 

The six baselines of the near detector(ND) and the far detector(FD) of 
RENO are L^ear (meters) = {660,445,302,340,520,746} and Lfar(meters) 
{1560,1460,1400,1380,1410,1480}, while their flux-weighted averages Inear and Ifar 
are 407.3m and 1443m, respectively. The baselines of Daya Bay, named EHl, EII2, and 
EH3, have lengths of EHl=494m, EH2=554m, and EH3=1628m, so that, conventionally, 
EHl and EH2 are regarded as near detectors while EH3 is regarded as a far detector. Since 
the energy of the reactor neutrinos falls mostly in the range of 2 MeV to 8 MeV, the typical 
L/E for reactor neutrino oscillation experiments like RENO and Daya Bay is estimated to 
be between 60 km/MeV and 270 km/MeV for the ND, and between 180 km/MeV and 800 
km/MeV for the FD. The coverage of the detectors of RENO and Daya Bay is described in 



Figffl 
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FIG. 1: The 4-neutrino oscillating aspects of Pth in Eq.© according to the baseline-to-energy 
ratio, L/E^ for different orders of Am|^. The shaded areas ND and FD denote the distances to 
the near detector and the far detector divided by the energy range 2 to 8 MeV. 
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FIG. 2: Four-neutrino analysis for the observed to expected ratios at both ND and FD: The ranges 
of i?ncar and i?far in the insets released by RENO and Daya Bay are reinterpreted in terms of the 
availability in sin^ 2^13 and sin^ 2^14 in a four-neutrino oscillation. 

The rate-only analysis of neutrino oscillation takes the average over accessible energies 
of the neutrinos emerging from the reactors. The measured probability of survival is 

^ j P,{E)a,,,{E)<t>{E)dE 
^ ^ j a,.,{E)<P{E)dE ' 

where (Jtot(-E') is the total cross section of inverse beta decay (IBD), and (t){E) is the neutrino 
flux distribution from the reactor. For Reno, the survival probability Pg in Eq.([8]) at each 
detector is evaluated as P„ear = 0.0678P„i + 0.1493P„2 + 0.3419P„3 + 0.2701P„4 + 0.1150P„5 + 
0.0558P„6 and P^r = 0.1373P/1 + 0.1574P/2 + 0.1809P/3 + 0.1856P/4 + 0.1780P/5 + 0.1608P/6, 
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based on the relative distances. Each Pni or Pfi is given as Pth in Eq.([7]). The total cross 



[0 



section of IBD is given as [26| 

= 0.0952 (^^^^^^"j X 10~V, (9) 

where = Ey — {Mn — Mp). The flux distribution 4>{E) from the 4 isotopes 
^U235 p^239 u238 p^24i^ at the reactors is expressed by the following exponential of a 5th 
order polynomials of E^, 

4>m=exp(j2fiE'}j, (10) 

where /o = +4.57491 x 10, /i = -1.73774 x 10-\ = -9.10302 x IQ-'^Js = -1.67220 x 
10-5, /4 = +1.72704 X 10-^ and /s = -1.01048 x 10"^ are obtained by fltting the total flux 
of the four isotopes with the flssion ratio expected at the middle of the burn up period of 
]. Including the product of crtotiE^) and (p{Ey) in the integrand in Eq. 



the reactors 
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results in average probability curve (P) shown in the inset of Fig. d When sin^ 26113 = 0.113 
(as announced by Reno), = curve is consistent with the released result. 

On the other hand, the probability at Daya Bay is evaluated by three groups of six 
reactors, EH1(AD1 & AD2), EH2(AD3), and EH3(AD1, AD2, & AD3), which catch neu- 
trinos from three groups of reactors, Dl & D2, LI & L2, and L3 & L4. The survival 
probability at each detector from different reactors is evaluated as follows: P[EH1) = 
0.795Pai + 0.143Pbi + 0.062Pci, P{EH2) = 0.065Pa2 + 0.512Pb2 + 0.423Pc2, and P{EH3) = 
0.246Pa3+0.379Pb3+0.375Pc3, where P^i, Psi and Pa are the probabilities at the i-th group 
of detectors due to reactors, Dl & D2, LI & L2, and L3 & L4, respectively. The coefficients 
are determined according to the relative distances. Substituting {P{EH1) + P{EH2))/2 
into Ps in Eq.([H]) gives the survival probability at near detectors. Likewise, substituting 
P{EH3) into P5 gives the probabihty at the far detectors. Although the ffux distribution 
at Daya Bay should have been obtained independently, (t){Ey) for RENO in Eq. lITOl) was 
included in the integration and, fortunately, the 6*14 = plot appears consistent with the 
published curve for sin^ 2^13 = 0.089 from Daya Bay, as shown in the inset of Figj2l 

In rate-only analysis, the value of sin^ 26'i3 is determined by obtaining the curve of (P) 
in Eq.(IH]) which can match the measured ratios Ps of the observed to expected ffux at each 
detector. However, there is a technical difference in the deffnition of the expected ffux using 
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RENO's approach compared to Daya Bay's. The definition of R that each experiment used 
is discussed in details in the next section. 

The main figures in FigJ2] interpret the measurements of the far-to-near ratio in terms 
of a four-neutrino oscillation. Two shaded bands in each frame indicate the ratios of the 
measured to expected at ND and FD, where the expected fiux is an estimation of the fiux 
expected without neutrino oscillation. The two bands in the plane of sin^ 2^13 and sin^ 2^14 
for the RENO experiment have ratios of i?far = 0.920±0.017 at FD and i?near = 0.989±0.025 
at ND. The intersection indicates the combination of sin^ 26is and sin^ 26ii allowed by the 
error bars in the inset. Meanwhile, the two bands for the Daya Bay experiment showed ratios 
of /?far = 0.931 ± 0.017 at FD and Kear = 0.986 ± 0.007 at ND; likewise, the intersection 
indicates the allowed combination of the two angles by the error bars in the inset. Depending 
on the existence of the 4th neutrino and the magnitude of sin^ 26'i4, the value of sin^ 26'i3 
might have a different value from that previously determined, as shown in Fig. |2l Such a 
combination of sin^ 2^13 and sin^ 2^14 will be examined more carefully within spectral shape 
analysis. 



IV. SHAPE ANALYSIS 



A. RENO 



One of the RENO's first results was the ratio of the observed to the expected number 

n 

of antineutrinos in the far detector, R = 0.920 ± 0.017 (see Ref. [5|), where the observed is 
simply the number of events at FD. On the other hand, the expected number of events at 
FD can be obtained through several adjustments of the number of events at ND, i.e., 

_ [Observed at FD] 

[Expected at FD] ^ ^ 

[No. of events at FD] , ^ 

- L J (12) 



[No. of events at ND]* ' 

where the number of events at each detector is normalized. The normalization of the neutrino 
fiuxes at ND and FD requires an adjustment between the two individual detectors which 
includes corrections due to DAQ live time, detection efficiency, background rate, and the 
distance to each detector. The numbers of events at FD and ND in Eg. (1121) have already 
been normalized by these correction factors, and so we have i?far = 0.920 ± 0.017 and 
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FIG. 3: The curves of the ratio of the survival at FD to the survival at ND of IBD antineutrinos: Red 
represents the A?ti|]^ oscillation in the three-neutrino analysis, while the blue represent the superpo- 
sition of two types of oscillations: one for Amj^ = 0.00232eV2 and the other for Am|i = O.OlGeV^. 
The data and errors in each figure are reproductions of original data. For the RENO(Daya Bay) 
setup, the red curve is drawn by Eq. i^{Eq. (dH)) with sin^ 26^13 = 0.124(0.0936) and the 
red dashed uses sin^ 2^13 = 0.113(0.089). The blue curves representing the superposition of two 
oscillations are obtained by sin^ 2^14 = 0.05 and 0.1, respectively. 

-Rnear = 0.989 ± 0.025 as shown in Figj2l The normalization guarantees i? = 1 at the center 
of reactors. RENO gets rid of the oscillation effect at ND when evaluating the expected 
number of events at FD by dividing the denominator of Eq. (fT2|) by 0.989 which is taken 

from i?near- NoW, 

^ [No. of events at FD] 

~ [No. of events at ND] / 0.989 ' ^ 

In rate-only analysis, the ratio of the observed to the expected number of events at FD in 

Eq.( !TT|) is just the survival at FD, since the denominator in Eq. ( IT3|l is eliminated. Thus, 

R coincides with the i?far in Fig. [2J 

In spectral shape analysis, however, the denominator cannot be neglected, since the 

oscillation effect at ND differs depending on the neutrino energy. The data points in Figj3] 

are obtained by the definition of the ratio R given in Eq. ffTTl) and Eq. f|T2l) per 0.25MeV 

bin, as the energy varies from 1.8MeV to 12.8MeV. The theoretical curves overlaid over the 

data are also obtained parallel to the ratio in Eq.( IT3|) . They are described by 



where PthiL) is given in Eq.([7]). While the thick red curve is a typical Am|^-dominant 
oscillation at 614 = 0, the blue curves are examples of a superposition of a Amg^^ oscillation 
and a Am|]^ oscillation. The red curve describes the oscillation related to Eq. ( 1141) when 
sin^26'i3 = 0.113, as determined from the rate-only analysis. However, if the data are 
compared with the curve of the ratio in Eq. iHM at On = 0, the minimum of is obtained 
at sin^2^i3 = 0.124, which is slightly different from sin^ 26'i3 = 0.113, the result of the 
rate-only analysis. The Ax^ with respect to sin^ 26'i3 is drawn in Fig. |5t^c). 

Even when nonzero 6'i4 is considered, the accessible range of Aml^ to catch the 4th 
neutrinos at RENO or Daya Bay is very narrow, ( just above Am|^ = 0.00232eV^), since 
the baselines are optimized for a Amg^^ oscillation. The interpretation of the data points 
in Fig. [3] in terms of the combined oscillations of Am|^ and Amli can be expressed as in 
Fig. m As expected, the Amli above 0.03eV^ is difficult to analyze in the four-neutrino 
oscillation scheme. The three curves in the figure correspond to the three exclusion curves 
at Icr, 2cr and 3a CLs as the result of shape analysis, i.e., with 99.7% confidence, values of 
sin^26'i4 larger than 0.3 are excluded. Figure [5] describes the preference of the combination 
of sin^26'i3 and sin^26'i4 for a given value of Am4^. Each figure contains la, 2a and 3a 
standard deviations, and the values of Am\^ for figures (a) and (b) were chosen from a 
couple of crests. Figure (c) shows a slope of Ax^ with respect to sin^ 26^13 at sin^ 26*14 = 0, 
which can be compared with the results of the three-neutrino rate-only analysis. 



B. Daya Bay 

In many aspects, the first results released by Daya Bay and RENO are parallel. There is 
a slight difference in what is expressed by the ratio R and the data points in Fig. [3] between 
the two experiments. Daya Bay also includes the result that R = 0.944 ± 0.008 (the ratio 
of the observed to the expected number of antineutrinos assuming no oscillations at the far 
detector), which implies that 

_ [No. of events at FD] 

" [No. of events at ND] ' ^ 

where the numbers of events at both detectors have been normalized for corrections including 
DAQ live time, detection efficiency, background rate, and distance to each detector. In 
comparison with RENO's R from Eq. fll2l) . the oscillation factor at ND is not eliminated 
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FIG. 4: The la, 2a and 3a exclusion curves. Due to the hmits of the basehnes and the neutrino 
energy, Amli larger than 0.05eV^ is excluded from the analysis. The range below O.OOSeV^ is 
blocked for Aml^, since oscillations in that range cannot be detected with a km-order baseline. 
For both, broad ranges of sin^ 2^14 apparently remain unexcluded. 




FIG. 5: Four-neutrino analysis in the sin^ 2^13 and sin^ 2^14 plane for chosen values of Am|]^: (a) 
0.028eV^ and (b) 0.016eV^. The last figure shows the best fit of sin^ 26*13, 0.124, and the standard 
deviation without sin^ 2^14 . 
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from the number of events in order to give rise to the expected number at FD. In other 
words, the expected number at FD takes the number of events at ND directly, while the 
observed at FD is the number of events at FD. Relative to the flux emerging from the 
reactors, the events at ND and FD are i?near = 0.986 ± 0.007 and Rf^, = 0.931 ± 0.017, so 
that 

R = 3^ = 0.944 , (16) 

which was the result released from rate-only analysis. 

The data points in Fig. [3] also gives R in Eq. (fT5|l per 0.25MeV bin, as the energy varies 
from l.SMeV to 8.5MeV. The theoretical curves overlaid for comparison with the data are 
consistent with the ratio in Eq. (fT5|) . and can be described by 

-Pth(-^far) ^^^-j 
near I 

which can be compared with the plots of Eq. (IT^ for RENO. While the thick red curve is 
a typical Amg^^-dominant oscillation without the 4th neutrino, the blue curves are examples 
of the superposition of a Amli oscillation and a Amli oscillation. The red curve describes 
Eq.f lT7|) when sin^2^i3 = 0.089 as determined from the rate-only analysis in Ref.^|, and is 
actually a reproduction of the curve from Ref. {4|. However, if the data are compared with 
the curve of Eq.( [T7|) at Ou = 0, the minimum of is obtained at sin^ 2^13 = 0.094, which 
is slightly different from sin^ 2^13 = 0.089, the result of rate-only analysis. The Ax^ with 
respect to sin^ 2^13 is drawn in FigEt^c). 

When nonzero 614 is considered, the interpretation of the data points in Fig. [3] in terms of 
the combined oscillation of Amg^^ and Aml^ can be expressed as in FigJH As in RENO, Amli 
above 0.02eV^ is avoided in the four-neutrino oscillation schemes. The three exclusion curves 
corresponding to la, 2a and 3a CLs are given from the result of spectral shape analysis, 
i.e., with 99.7% confidence, values of sin^ 2614 larger than 0.3 are excluded. For a given value 
of Am4^, the combination of sin^ 2^^i3 and sin^2^i4 is analyzed in Fig. El where each figure 
contains la, 2a and 3a standard deviations. In figures (a) and (b), the values of Amli are 
chosen to be 0.032eV^ and 0.016eV^ for a couple of crests in FigJU respectively. Figure (c) 
shows the slope of Ax^ with respect to sin^ 26'i3 at sin^ 2^14 = 0, where Ax^ is minimized at 
sin2 20i3 = 0.0936. 
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FIG. 6: Four-neutrino analysis in the sin^ 2^13 and sin^ 2^14 plane for chosen values of Am|]^: (a) 
0.032eV2 and (b) 0.016eV^ The last figure shows the best fit of sin^ 26^13, 0.0936, and its standard 
deviation without sin^ 2^14. 

V. CONCLUSION 

If a fourth type of neutrino has a mass not much larger than the other three, the results 
of reactor neutrino oscillations like RENO, Daya Bay, and Double Chooz can be affected by 
the fourth state. For detectors established for oscillations driven by Amg^ = 0.00232eV^, 
clues of the fourth neutrino can be perceived only if the order of Amli is not much larger 
than that of Amg^^. Therefore, this work examined the possibility to find a kind of sterile 
neutrino for the range of mass-squared difference below O.OSeV^. On the other hand, value 
of Amli below O.OOSeV^ is not considered either, so that an approximation was used in 
evaluating oscillation probabilities. Otherwise, the contribution of Am|j^ oscillation to the 
disappearance of IBD neutrinos is difficult to separate from that of Am|]^ oscillation. We 
examined the two announced results of RENO and Daya Bay, in terms of a four-neutrino 
oscillation for a certain range of Am^^^. 

The first results released by the two experiments are rate-only analyses providing the 
far-to-near ratio of properly normalized events. The ratios, R = 0.920 ± 0.017 from RENO 
and R = 0.944 ± 0.008 from Daya Bay were interpreted as sin^ 26113 = 0.113 ± 0.023 and 
sin^ 2^13 = 0.089 ±0.011, respectively, in their original three-neutrino analysis. For compari- 
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son, we interpret the ratios with respect to a four-neutrino oscillation. The range of sin^ 2^13 
broadens along with the range of sin^ 2614, as shown in Fig. |2l which took (and shifted) the 
errors from the i?'s. 

Although the spectral shape analysis was not presented due to lack of data in the first 
releases of RENO and Daya Bay, we included the shape analysis of each oscillation as 
shown in FigJS] Because of the accessibility of the baseline, Aml^ larger than O.OSeV^ is 
excluded. As expected, only Am|^ below O.OSeV^ exhibits modest oscillatory aspects for 
both in FigJH Regarding RENO, Fig. [3], which examines the exclusion boundary in the 
Am^i-sin^ 2^14 plane, Aml^ = 0.042eV^ and 0.028eV^ seem the most likely for arbitrary 
sin2 2^i4. When Amj^ = 0.042eV^(0.028eV^), values of sin^ 20i4 larger than 0.43(0.53) are 
excluded by Sex CL as shown in FiglSl Regarding Daya Bay, there is no preferred Amli 
and values of sin^ 26'i4 > 0.38 are excluded by 3cr CL as shown in FigJH When sin^ 26'i4 is 
taken into consideration, the minimum is given rise to by sin2 2^i3 = 0.124 for RENO 
and sin^ 2^13 = 0.0936 for Daya Bay, which are slightly different from the rate-only analysis 
results of the first announcements. In this rough estimation, sin^ 2^13 = is not excluded 
by 3a CL for either RENO or Daya Bay. 
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